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Abstract 

Radiation reaction (but, more generally, fluctuations and dissipation) occurs when a system in- 
teracts with a heat bath, a particular case being the interaction of an electron with the radiation 
field. We have developed a general theory for the case of a quantum particle in a general poten- 
tial (but, in more detail, an oscillator potential) coupled to an arbitrary heat bath at arbitrary 
temperature, and in an external time-dependent c-number field. The results may be applied to a 
large variety of problems in physics but we concentrate by showing in detail the application to the 
blackbody radiation heat bath, giving an exact result for radiation reaction problem which has no 
unsatisfactory features such as the runaway solutions associated with the Abraham-Lorentz theory. 
In addition, we show how atomic energy and free energy shifts due to temperature may be calcu- 
lated. Finally, we give a brief review of applications to Josephson junctions, quantum statistical 
mechanics, mesoscopic physics, quantum information, noise in gravitational wave detectors, Unruh 
radiation and the violation of the quantum regression theorem. 
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I. INTRODUCTION 



Radiation reaction is familiar to most readers through the Abraham-Lorentz equation for 
a radiating electron [1-3] . It arises from the fact that electric and magnetic fields emitted by 
an accelerating electron act back on the electron, resulting in a retarding force. However, the 
equation obtained by Abraham [1] and Lorentz [2] displays unphysical runaway solutions. 
Through the years, there have been many attempts to obtain a more satisfactory result but 
none were without difficulties of one kind or another until in 1991 we proposed a physically 
consistent solution in the form of a quantum Langevin equation [4] . We discuss that result 
below in Sec. Ill, but first we consider more general properties of dissipation and fluctuations, 
since they appear in many areas of physics. 

The equation of motion of a quantum particle in a heat bath (also referred to as a reservoir 
or as the environment) was pioneered in 1965 by Ford, Kac and Mazur [5]. These authors 
used a quantum model of coupled oscillators to obtain a quantum Langevin equation that 
showed explicitly the role of fluctuation and radiation reaction (dissipation) effects. By 
1983, there was much general interest in mesoscopic systems and Caldeira and Leggett [6] 
showed that quantum heat baths could play an important role in the analysis of such systems 
by analyzing in detail the effect of dissipation on quantum tunneling in such systems as a 
single Josephson junction or a SQUID. These authors, again considering an Ohmic heat bath, 
used the path integral influence functional approach of Feynman and Vernon [7], a technique 
which was then followed by many investigators, particularly these working on coherence and 
entanglement problems. By contrast, the present authors, in collaboration with J. T. Lewis, 
developed a general quantum Langevin equation approach to these problems [8, 9] that, 
in our opinion, is simpler and more physically transparent. In particular, we went beyond 
the Ohmic model to consider an arbitrary frequency- dependent heat bath (in particular a 
blackbody radiation field heat bath). 

In general, coupling a system to a heat bath produces two related effects: dissipation 
and fluctuations. The study of these phenomena is loosely referred to as the "Brownian 
motion," problem, going back to the work of Robert Brown, a Scottish botanist who in 1828 
and 1829 published his work on the random motion of pollen grains immersed in a fluid 
[10]. The explanation of Brown's phenomenon was given by Einstein and Smoluchowski in 
the first decade of the twentieth century [11, 12]. This work helped to definitively establish 
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the atomic theory of matter since the irregular motion was clearly identified as being due 
to collisions with the molecules in the liquid. The term "Brownian motion," is now used 
in a generic sense to denote random motion and it covers a wide spectrum of phenomena 
from the motion of very fine particles suspended in a gas to the motion of electrons in a 
blackbody radiation heat bath. 

Einstein's explanation of Brownian motion used a discrete time approach. In particular, 
his results included an explanation of the fact that increasing temperatures lead to more 
agitated Brownian motion, in a relation in which the diffusion coefficient is shown to be 
equal to the product of kT and the mobility, perhaps the first example of what is now called 
the fluctuation-dissipation theorem. 

Shortly after the work of Einstein and Smoluchowski, Langevin [13] presented a different 
approach to the subject which, in the words of Chandrasekhar [14], constitutes the " mod- 
ern," approach to this and other such problems. The essence of Langevin's approach is a 
continuous time approach implemented by adoption of a stochastic differential equation, i.e., 
an equation for quantities which are random in nature. In other words, Langevin provided 
an elegant solution to the problem of generalizing a dynamical equation to a probabilistic 
equation. This was to be the start of a major new field of study with applications in physics, 
chemistry, biology, and many other fields. 

The approach of Langevin was phenomenological but its essential correctness has been 
verified by various microscopic studies. A key feature of his approach was to separate the 
total force acting on a particle due to its environment into two parts: a mean frictional force 
and a fluctuation (random) force. These forces are related to each other as a consequence of 
the requirement that the equilibrium state be stationary. On the other hand these two forces 
are very different in nature: The fluctuation force is basically microscopic in nature and has 
a time scale determined by the mean time between collisions, whereas the time scale of the 
frictional force is proportional to the self-diffusion constant and is much larger. Explicitly, 
Langevin explained the motion of a " Brownian," particle (an otherwise free particle in a 
dissipative environment) by an elegant stochastic classical differential equation [13]: 

mx + (x = F(t), (1) 

where m and x denote the mass and coordinate of the particle, respectively, and the dot 
denotes differentiation with respect to time. The force on the particle consists of the frictional 
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(dissipative) term — (x and the random (fluctuation or noise) term F(t). The latter term is 
zero at zero temperature, in contrast to the corresponding situation in quantum mechanics, 
as we shall see later. 

Since the past motion does not appear in the Langevin equation (1), one says that 
there is no memory or, equivalently (in the classical case under discussion but not in the 
quantum case to be discussed later), we say that the process is Markovian. In addition, 
the autocorrelation of the random force is a 5 function (a manifestation of the Markovian 
process) and is also proportional to (. The latter result is a manifestation of the fluctuation- 
dissipation theorem. 

In addition, with the usual definition of the diffusion constant, 



is the mean-square displacement and the angular brackets denote the average with respect 
to the canonical ensemble of the system, Langevin readily obtained the famous Einstein 
relation 

*-f. 

As noted above, this is an example of an intimate connection between fluctuations and 
dissipation. Another example was provided by Nyquist who showed that the random fluc- 
tuations in voltage across a resistor measured by Johnson are determined by its impedance 
(the famous Johnson-Nyquist noise [15] in electrical circuits). A general quantum formula- 
tion of the fluctuation-dissipation theorem first appeared in the well known article of Callen 
and Welton [16]. 

We turn now to another key property of the Langevin equation: Whereas the original 
Brownian motion experiment and the analysis of it pertained to a system in thermal equi- 
librium, the Langevin equation can be generalized in a simple manner to include an external 
force on the right hand side and thus can describe the irreversible approach to equilibrium. 

The Langevin equation is phenomenological. As noted above, its derivation from a quan- 
tum microscopic theory was first given by Ford, Kac and Mazur [5]. Then in 1986, in an 
article entitled "Quantum Langevin Equation", [9] together with J. T. Lewis we gave a de- 
tailed discussion and presented the general form of the equation consistent with fundamental 




(2) 



where 



s(t-t>) = ([x(t)-x(t>)} 2 ) 



(3) 
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physical requirements, in particular causality and the second law of thermodynamics. Based 
on these conclusions, we were able to obtain the most general quantum Langevin equation 
[given below in (6)] for the macroscopic description of a quantum particle with passive dissi- 
pation and moving in an arbitrary external potential. We then showed that the most general 
form can be realized by a simple oscillator model of a heat bath. For the purpose of this 
article, it is convenient to use this model, the so-called independent-oscillator (10) model, 
for which the quantum Hamiltonian is 

Hio = ^ + V(x) + £ (JL + \m jU * {q 3 - xf\ - xf (t) . (5) 

j \ 3 / 

Here, V(x) is a particle potential energy, / (t) is an external applied force, while x and p 
are the particle position and momentum operators and q 3 - and p 3 are those for the j'th bath 
oscillator. The parameters are m, the particle mass, and m 3 and Uj, the mass and frequency 
of the j'th bath oscillator. 

The procedure used in obtaining the corresponding Langevin equation from the micro- 
scopic Hamiltonian is common to all such problems and consists of use of the Heisenberg 
equations of motion to obtain the equations of motion for both the dynamical variables 
of the particle (x,p) and the dynamical variables of the heat bath (qj,Pj) [9]. These are 
coupled equations and the next step is to eliminate the bath variables. This leads to an 
inhomogeneous differential equation for the q 3 . Then, typical of the way that the time- 
reversal invariance of the original equations is broken in macroscopic equations, one chooses 
the retarded solution of this equation. This solution for q 3 is then substituted into the equa- 
tion for x to get the Langevin equation for a quantum particle of mass m moving in a one 
dimensional potential V(x) in an arbitrary heat bath and temperature T: 

mx+ j dt 1 n(t-t x )x(t x ) + V'(x) = F{t) + f{t), (6) 

J — oo 

where the dot and prime denote, respectively, the derivative with respect to t and x. This 
is a Heisenberg equation of motion for the coordinate operator x. The coupling with the 
heat bath is described by two terms: an operator-valued random force F(t) with mean zero, 
and a mean force characterized by the memory function fj,(t). These quantities are given in 
terms of the heat bath variables: 

H(t) = ^rrijUj] cosiujj^eit), (7) 
j 
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where 6(t) is the Heaviside step function (by convention the memory function vanishes for 
negative times), and 

F(t) = J2m J ^(t), (8) 

3 

where qj(t) denotes the general solution of the homogeneous equation for the heat bath os- 
cillators (corresponding to no interaction). Using these results, we find that the (symmetric) 
autocorrelation of F(t) is 

- (F(t)F(t f ) + F(t')F(t)) = - f duRe + i0 + )] hu coth cos[w(t - *')], (9) 
2 7r Jo 2/cT 

and the nonequal-time commutator of F(t) is 

[F(t), F(t')} = — / duRe \fi(u + i0 + )} hu sm[u(t - if)]. (10) 
In these expressions 

POD 

jl(z) = / dte izt fi(t), Imz > 0, (11) 
Jo 

is the Fourier transform of the memory function /i(t). Finally, F(t) has the Gaussian prop- 
erty: correlations of an odd number of factors of F vanish; those of an even number of factors 
are equal to the sum of products of pair correlations, the sum being over all pairings with 
the order of the factors preserved within each pair. Equation (9) is an exact fluctuation- 
dissipation theorem and we emphasize that it is independent of both the potential V(x) and 
the external force fit). Physically, it expresses the fact that the spontaneous equilibrium 
fluctuations of the heat bath (described by the left side of (9)) are related to the dissipation 
parameter Rejl(uj). 

From the explicit expression (7) it is clear that the memory function is dependent only on 
the bath parameters. We obtained a key constraint on this function by considering the effect 
of an arbitrary c -number external force f(t) acting on an otherwise free particle (V (x) = 0) 
[9]. If we assume that / (t) vanishes in the distant past and future, the effect is to carry the 
system of free particle coupled to the bath in a complete cycle from a state of equilibrium, 
through a continuous sequence of intermediate states, and back to equilibrium. The second 
law of thermodynamics, in the Kelvin-Planck form, then requires that the net work done by 
this force be positive, which in turn requires that the spectral distribution Re [fl [oj + i0 + )] 
must satisfy the positivity condition: 

Re{/2(w + i0 + )} > 0, -oo < u < oo. (12) 
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We note that since the memory function is real the spectral distribution must also satisfy 
the reality condition: Re [fx (—u + i0 + )] = Re [fx (u> + i0 + )]- 

The positivity condition (12) together with the fact, obvious from the definition (11), 
that fx(z) is analytic in the upper half plane, means that means that fx(z) must be what 
is called a positive real function. This is a very restricted class of functions of a complex 
variable, with special properties which include the Stieltjes inversion theorem: 

fa) = — rjj&p+pi. (13) 

vr J z 2 - u 2 

Thus, we see that in the quantum Langevin equation the memory function as well as the 
correlation and commutator of the random force, are completely characterized by the spectral 
distribution. 

There are further conditions on the spectral function. As a consequence of the inversion 
theorem (13) we see that 

f°° , Re{fi{uj + tO + )} 

I d " iW <oc - (14) 

An important further constraint is what has been called the zero'th law of thermodynamics: 
there must be an equilibrium state. This requires that the spectral distribution must be 
everywhere positive, with no gaps in which it vanishes. This in turn requires for our model 
that the bath frequencies must be infinite in number and continuously distributed. 

In the case of an oscillator potential V(x) = t\Kx 2 , the solution of the quantum Langevin 
equation (6) is given by [8] 



x(u) = a(Lu 



F(u) + f(u) , (15) 



where the superposed tilde to denotes the Fourier transform, e.g., x(uj) is the Fourier trans- 
form of the operator x(t), and a(u) is the generalized susceptibility (a c-number) given 
by 

= 2 • ~t \ \ is - ( 16 ) 

—moj 2 — ioufj,(uj) + K 

As a simple application of this result, we find for the autocorrelation and commutator of 
the operator x(t), in linear response, are [8, 17] 

1 Ti f°° hco 

- (x(t)x(t') +x(t')x(t)) = - dculm{a(u + iO + )}coth—— cos (to (t-f)), 

2 7r J 2kT 

[x(t),x(t')} = — / dcolm{a(uj + i0 + )} sin u (t - t'). (17) 
m Jo 



Equations (9), (10) and (17) are exact fluctuation-dissipation theorems for our general 
analysis and they provide the foundation for subsequent developments. We note, however, 
that equation (17) depends on both the potential and the dissipative parameter. In fact, 
in the classical case (but not in the quantum case), we showed that, in the case of a weak 
applied force, the spontaneous fluctuations (described by the left-side of (17)) relax with the 
same time constant as the induced (by fit)) non-equilibrium fluctuations [18]. This is the 
essence of the Onsager regression hypothesis [19] which states that regression of fluctuations 
is governed by macroscopic equations describing the approach to equilibrium. 

We have now essentially all the tools that we need and thus it is time to turn to specific 
applications. However, here we decided to concentrate on important problems in electrody- 
namics. First, in section II, we survey the vast amount of work on the problem of runaway 
solutions of the Abraham-Lorentz equation for a radiating electron. Many of these attempts 
started with force equations which do not encapsulate the important time development na- 
ture of the problem. By contrast, our starting point is a Hamiltonian which beautifully leads 
to an equation of motion valid for all times. Moreover, they treated the acceleration x(t) as 
the key parameter, ignoring the fact that it is not an observable, in contrast to the applied 
external force f(t) (which, of course, differs from Mx(t) because of the dissipation), which 
is the basic parameter used in our approach. Since all previous solutions suffered from one 
problem or another, in section III we present our consistent theory, which led us to a simple 
second-order differential equation which is not only free of runaway solutions but has no 
causality problems and is consistent with the optical theorem. Since the basic framework 
is the same, we also include in this section our treatment of temperature effects on atomic 
energy and free energy levels. 

In section IV, we present the relativistic extension of our theory and in section V we 
discuss associated fluctuation and quantum effects. Then, in section VI, we briefly enumer- 
ate other applications of our general theory which embraces both radiation reaction and 
dissipation. Our conclusions are summarized in section VII. 
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II. RADIATION REACTION IN ELECTRODYNAMICS: HISTORICAL SURVEY 



The earliest work on radiation reaction is that of Abraham [1] and Lorentz [2], whose 
result is summarized in the well-known equation: 

Mx - Mr e x = f(t) Abraham-Lorentz, (18) 

where 

-e = - 6.25 x 10-» 8 . (19) 

For a careful discussion of the derivation of this equation, which is "exact," for a point elec- 
tron, see Jackson [3]. This equation exhibits the well-known problem of runaway solutions: 
even a small impulsive force acting on an electron at rest results in an exponentially growing 
displacement. This is made more explicit in the sentence following equation (45). In fact, 
when f(t) = 0, equation (18) does not reduce to Newton's equation, as it should. In essence, 
the problem with the equation might be thought to lie in its derivation which was based on 
force equations, as distinct from a Hamiltonian. However, as we shall later point out, our 
exact Hamiltonian approach leads, in the case of a point electron to the same equation (18), 
making clear that the basic assumption of a point electron is at fault. In fact, both Abra- 
ham and Lorentz presented a more systematic discussion by considering both the charge 
structure of the particle and its self-fields. However, their derivation contained a variety of 
assumptions and their final result was an infinite expansion, the leading term corresponding 
to the point electron result given in (18). 

Attempts to solve this problem have engaged the efforts of a large number of investigators 
over the past century. Since, in our view, none of these efforts have been successful, we will 
just concentrate on presenting the more prominent. First, we note that Born and Infeld [20] 
attempted to fix the problems associated with the Abraham-Lorentz approach by modifying 
Maxwell's theory to make it non-linear but, in particular, they encountered problems with 
quantization. 

Dirac [21] attempted to solve the problem by including advanced solutions in addition to 
retarded solutions but this effort suffered from a violation of causality (so that the accelera- 
tion at time t depended on the force acting at times earlier that t). An attempt by Ivanenko 
and Sokolov [22] to replace the Abraham-Lorentz equation by an integro-differential equa- 
tion of motion was also flawed by virtue of introducing a violation of causality. Similar 
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remarks apply to the theory of Wheeler and Feynman [23]. See also [24-26]; we discuss 
the attempts of Eliezer [25] and Landau and Lifshitz [26] in more detail in section III, in 
connection with our own results. 

III. A CONSISTENT THEORY OF RADIATION REACTION 

The Abraham-Lorentz equation is a result of nineteenth century physics with, moreover, 
no notion of fluctuations. Almost without exception, discussions of radiation reaction in the 
past century have relied on the same physics. On the other hand, our contribution has been 
to recognize that a correct equation is a quantum Langevin equation for which the heat 
bath is the blackbody radiation field, with fluctuations due to the fluctuations of that field. 
As we have seen, a consistent derivation of such an equation begins with a Hamiltonian 
formulation of the dynamics and the laws of thermodynamics impose a powerful constraint 
on that formulation. 

For a nonrelativistic electron (charge = — e) interacting with the quantum electrodynamic 
radiation field the Hamiltonian has the form [9]: 

Hqeb = ^ (p + ^A) 2 + V(r) + ^4,s«M - xf (f) , (20) 
where the vector potential is given by 




Here the symbols have their usual meanings. The quantity fk is the electron form factor 
(Fourier transform of the electron charge distribution). Without loss of generality, we have 
taken the form factor as well as the polarization vector e^, to be real. The form factor, 
which is sometimes called a cutoff factor, must have the property that it is unity up to some 
large cutoff frequency f2 beyond which it falls to zero. Then, we showed [9] that, by a unitary 
transformation that leaves the position operator unchanged, this Hamiltonian can be put 
into the independent oscillator form (5). Rather than repeating that discussion, we will show 
how to get the quantum Langevin equation directly from this form of the Hamiltonian. 
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The Heisenberg equations of motion are 

in m 
1 dV 
P = tt/[p,#qed] = — + 



1 / 27re 2 p + -A 

m y TiUkV m 

Eliminating the particle momentum operator between the first two equations, we get the 
particle equation of motion: 

mr+^- = -A. (23) 

OT C 

The solution of the last equation can be written 



«M (t) = a£> (f) - i\^^h^,s ■ f ^ dt'e-^-^r (?) , (24) 

where aj^] (t) is the solution of the homogeneous equation, corresponding to free motion of 
the bath in the absence of the electron. Putting this in the expression (21) for the vector 
potential, we find 

A (t) = A^ (t) - ^ — f dt ' sin N (t ~ t') e Ks e Ks • r (*')] . (25) 

The sum over s is the sum over the two polarization directions perpendicular to k, so we 
have Yls &k,s&Ks • r = r — k • rk. Next we form the limit of an infinite quantization volume, 
using the prescription — > -^-^ J dk. With this the particle equation of motion (23) takes 
the form of a quantum Langevin equation: 

r* dv 

mr + J dt'n (t - t') r (*') + — — F(t) + f (t) , (26) 
where the memory function is 

4e 2 f°° 

V (t) = ^ J du k u 2 Jl cos (ou k t) 9 (t) , (27) 
and the fluctuating operator force is 

F (t) = -k {h) it) . (28) 

c 

Note that these expressions for the memory function and the fluctuating force satisfy the 
general feature that they depend only on the bath parameters, independent of the particle 
mass and the external potential. 
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From the expression (27) for the memory function, we find that the spectral distribution 

is 

Re [H (u + i0+)] = ^ J du k uif k 5 (co - u k ) = 2 -^f 2 . (29) 

The physically significant results for this model should not depend upon details of the 
electron form factor, subject, of course, to the condition that it be unity up to some large 
frequency Q and falls to zero thereafter. A convenient form which satisfies this condition is 

* - z&r < 3 °> 

Using this in the expression (29) for the spectral distribution, the Stieltjes inversion formula 
(13) gives 

-/ n 2e 2 ztt 2 

M = (31) 

Finally, with the form (30) for the form factor, the expression (27) for the memory function 
can be evaluated to give 

fi (t) = ^n 2 (25 (t) - ne' nt ) e (t) . (32) 

Put this in the equation of motion (note S(t)6(t) is "half S(t)) with 
V (r, t) = \Kr 2 — r • f (t) , corresponding to an external oscillator potential and an applied 
force f {t). Then multiply both sides by e nt and differentiate with respect to t, to get the 
equation of motion in the form: 

tfnr + (m + r + ±tfr + Kr = f (t) + if (t) +F(t) + (t) . (33) 

For motion for which the typical frequencies are much below the cutoff frequency Q, this 
becomes of the form for a free particle, but with a renormalized mass: 

2e 2 Q . . 

The mass M is interpreted as the observed mass of the electron, although at ultrahigh 
frequencies m, the bare mass, reappears. Now with this interpretation the bare mass, what 
we have up till now been calling the electron mass, is no longer observable, but is given in 
terms of the observed quantities through 

m = M (1 - ttr e ) , (35) 
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where M — 9.11 x i(T 28 # and r e = 2e 2 /3Mc 3 = 6.25 x l(T 24 s. With this expression for the 
bare mass, the equation (33) can be written 

M (j}~ ^ r + Mr+^Kr + Kr = F (t) + (t) + f (t) + if (t) . (36) 

Within the framework of our model, this is an exact Heisenberg operator equation of motion. 

For the purpose of making contact with the Abraham-Lorentz and other classical equa- 
tions, in Eq. (36) we will take mean values, set K — and specialize to one dimension. 
Thus, the fluctuating force F(t) is eliminated and all quantities are now to be interpreted 
as classical quantities and we obtain 

M (IT 1 - r e ) x(t) + Mx(t) = fit) + VI' 1 j\t). (37) 

We immediately see that the Abraham-Lorentz equation (18) follows if we take Q — > oo, 
corresponding to a point-electron model for the electron. But in this limit the bare mass 

(35) is negative infinity. This is the source of the runaway solutions that plague that model. 
If the bare mass is to be positive the relation (35) puts a constraint on the cutoff frequency: 

n < r- 1 (38) 

and hence the point electron model, associated with the Abraham-Lorentz equation, is 
ruled out. The largest possible value of the cutoff consistent with this constraint is = r" 1 , 
corresponding to zero bare mass. Choosing this value of the cutoff, the equation of motion 

(36) becomes 

Mr + r e Ki + Kr = f(t) + tJ (t) + F (t) + r e F (t) . (39) 

This is a rather striking exact result in that it is only a second-order equation with the only 
parameter being r e . Its form is a result of our choice (30) for the cutoff function with the 
choice Q = r" 1 . Other forms of the cutoff function will give rise to terms on the right hand 
side of higher order in r e , but the first order term is the same for all. Since these higher order 
terms reflect meaningless details of the cutoff function, we feel that the simple equation (39) 
is the one of choice. It is our key result. 

The case usually discussed in the literature corresponds to a free classical particle (K = 0) 
with neglect of fluctuations (F (i) = 0) and motion in one dimension, where our equation 
(39) specializes to the form 

Mx = f (t) + T J (t) . Ford-O'Connell. (40) 
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We emphasize that f(t) is a general time-dependent external field. In the particular case 
where this field is an electric field, we found that [27, 28] 

'MirWr e£(t) ' (41) 

which, except for very large frequencies u is essentially e E{t). It is of interest to note that 
Eliezer [25], whose approach was to postulate what he considered possible solutions, also 
wrote down (40) with f(t) = e E(t) [see Eliezer's equation (9)]. Landau and Lifshitz [26] 
later obtained by a method of successive approximations to the Abraham-Lorentz equation 
a result similar to that of Eliezer. Of course, neither of these investigators realized that all 
departures from the AL equation require the existence of a charge with structure. 

Equation (40) is to be compared with the Abraham-Lorentz equation (18). Indeed, if 
in that equation we assume that in first approximation the term —Mr e x is small and can 
be neglected, we get Mx = f (t) and —Mr e 'x = —r e f(t), giving our equation (40). This 
is essentially the argument used by Eliezer [25] and by Landau and Lifshitz [26] to get a 
corresponding result. But, due to the existence of runaway solutions, the term —Mr e 'x is 
never small for all times, so their argument is flawed. On the other hand, as we have seen, 
our equation (40) is the result of a consistent theory with consistent approximations. 

Forming the Fourier transform of equation (39), we can write the solution as 

r(u>) = a (u>) [f (u) + F (cu)] , (42) 
where the polarizability is given by 

aM= -M^ + (l-^r e )ir - (43) 
As a simple application we can calculate the autocorrelation (17) in the classical limit 
(k — > 0). With the polarizability (43) we find 

1 2kTr f°° to 2 

- (x(t)x(0) + x(0)x(t)) = -ri- 1 / du ^ cos (cut) 



' K 



1 

cos {u\t) sin (oj\ \t\ 

1oj\ 



(44) 



where 

u,l = . 

M ' M 



2 K Kr e / 9 7 2 / . _ 

UJ o = T7^ 7 = nrr> ^ = \ u o--r- (45) 
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If we make the same calculation with the Abraham-Lorentz equation we find an autocorre- 
lation that grows exponentially in time, a clearly unphysical result which shows that in the 
presence of fluctuations the problem of runaway solutions is inescapable. 

It is well known that the Abraham-Lorentz equation is compatible with the Larmor 
formula for the radiated power. However, for our equation (40) we find that the total 
electromagnetic energy radiated from a confined current distribution [28, 29] is 



This is our generalization of the Larmor formula. If we make the replacement f (t) — > 
Mr, which corresponds to setting r e = Q' 1 — > 0, we get the familiar form of the Larmor 
formula. Thus we see that our expression (46) for the radiated energy is compatible with 
our form (40) of the equation of motion, just as the Larmor expression is compatible with 
the Abraham-Lorentz equation. Finally, we should emphasize that for a force f (t) that is 
slowly varying on a time scale r e the difference between the two expressions is negligibly 
small but our expression is in terms of the aplied force, as distinct from the theoretically 
derived acceleration. 

At first, one might be surprised that (39) predicts that for a free particle (K = 0) there 
is no radiation if the external force is constant whereas (46) seems to state otherwise. The 
answer is that (46) was derived, as is usual, with the assumption that f(t) is switched on 
in the distant past and off in the distant future. For a force that is switched on, is constant 
for a long time, and then switched off, there is no radiation during the intermediate times of 
constant force [29], yet the total radiated energy is correctly given by the formula (46). This 
was verified in an explicit example [30]. With some exceptions, it appears that this simple 
explanation was often missed in the endless debate surrounding this problem in the past 
(where, in dealing with the Abraham-Lorentz equation, the argument centred on constant 
acceleration). The most notable exception was Feynman [31] who states that "- - we have 
inherited a prejudice that an accelerating charge should radiate - the power radiated by an 
accelerating charge [the Larmor formula] has led us astray," and he then goes on to discuss 
the limited validity of the Larmor formula and the fact that "- - it does not suffice to tell us 
" when," the energy is radiated" . 

Finally, we point out that our equation (39) is not only free of runaway solutions but 
is also consistent with the optical theorem and the standard formulas for the Rayleigh and 




(46) 
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Thomson scattering cross sections [27]. Moreover, the corresponding polarizability (43), 
since it is analytic in the upper half u plane, is consistent with the basic physical requirement 
of causality. 

It is also interesting to note that in our first paper in this general area [8], we considered 
the motion of a charged particle in a radiation field, with the purpose of finding the effect of 
temperature on atomic energy levels (and, of course, this was the framework we later used 
for the treatment of radiation from the electron). This followed much earlier work by Knight 
[32] at a time when it was considered that the effects were so small as to be unobservable. 
However, using high-precision laser spectroscopic techniques, in a remarkable paper, Hollberg 
and Hall [33] were able to measure fractional shifts of ~ 2 x 10~ 12 . In our approach, since 
temperature was involved, we were obliged to consider thermodynamics and, in particular, 
free energy. As a result, we obtained a striking exact result for the free energy of a quantum 
oscillator interacting with a blackbody radiation field, which we used to obtain agreement 
[34] with the experimental results obtained for the energy shifts due to temperature [33]. 
Later, when a flurry of paper appeared claiming the laws of thermodynamics were invalid 
in the quantum arena, we were able to show, using this same free energy result, that these 
claims were in fact incorrect [35]. 



IV. RADIATION REACTION: RELATIVISTIC THEORY 

Analogous to Dirac's extension of the Abraham-Lorentz equation [21], we have proposed 
an extension of our equation (39) to the relativistic domain [30]. The form is 

M *£ - l F >* + {l F > x - h^'l F ^) • < 47 > 

where F^" is the external electromagnetic field tensor and 



i / v 2 

u" = —, dr = -y%^dx^ = Wl - — dt. (48) 

This equation is consistent with the constraint 

g^uW = c 2 , (49) 

and in the nonrelativistic limit (c — > oo) reduces to our equation (39) for the special case of 
a free particle (K = 0) and for the mean motion with no fluctuation (F (t) = 0). Also in 
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that limit, the applied force is f (t) = eE(t) , which tells us that the equation is really only 
valid for spatially uniform (although possibly time dependent) fields. 

It is also of interest to note that our equation of motion (47) can be written in the 
three-vector form 



is the Lorentz force with E and B the electric and magnetic fields. We note that the 
corrections to the non-relativistic results are of order (f/c) 2 , as one might expect. 

In [30] we presented an exact solution of this equation for the case of an electron travelling 
between the plates of a parallel-plate capacitor, for which the electric field is uniform between 
the plates and zero otherwise. The result shows explicitly the radiation occurs only as the 
electron enters and leaves the field. 

V. FLUCTUATION AND QUANTUM EFFECTS 

Earlier discussions of radiation reaction were generally based on classical electrodynam- 
ics and implicitly assumed that fluctuations could be neglected. For the most part, our 
discussion above was equivalent with such a classical description, with the fluctuation force 
eliminated by taking mean values. As we have seen in our example (44), in the classical 
limit this corresponds to zero temperature, which is not a serious limitation for the classical 
theory. However, in the quantum theory zero point fluctuations are always present. As 
an example of the importance of quantum fluctuations, we consider the commutator of the 
position and velocity. Forming the commutator of r with the first of equations (22), using 
the canonical commutation we find when specialized to one dimension, [x(t),x(t)\ = ih/m, 
where m is the bare mass. However, using the expression (43) for the polarizability, the 
formula (17) for the non-equal time commutator can be readily evaluated and it is obvious 
that only the renormalized mass appears. In particular, in the limit as t' approaches t from 
above or below we find that [x(t),x(t±0 + )] = i-jfe (1 — 7r e ), with M the renormalized mass. 
For a detailed discussion see [36]. 




(50) 



Here v = r, 7 = \J\ — v 2 /c 2 and 




(51) 
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As another example of a quantum fluctuation phenomenon, consider the mean square 
displacement of a harmonically bound electron at zero temperature, obtained by setting 
t' — t in the correlation (17) and forming the zero temperature limit. With the expression 
(43) for a (oj), we obtain an expression in terms of an integral that is logarithmically diver- 
gent. However, inclusion of retardation would make the integral finite with an upper limit 
approximately equal to Mc 2 /%. [37] With this mean square displacement is found to be 



Since ojQT e is presumed very small this corresponds to a small increase over the leading factor, 
which is the mean square width of the oscillator ground state. 

Finally, we remark on quantum tunneling in a dissipative system. For most systems of 
interest it was found that dissipation decreases tunneling rates. However, in the case of the 
blackbody radiation field, we found that tunneling actually increased [38]. The reason for 
this exception to the general rule is the presence of mass renormalization. 

VI. MISCELLANEOUS APPLICATIONS 

The case of constant friction is of special interest. There the spectral distribution is 
independent of oj: Re [p, (oj + i0 + )] = (, the friction constant. This is frequently referred to 
as an Ohmic heat bath. The equation (6) then takes the form: 



This is the same form as the original, classical form of the Langevin equation but here, of 
course, x and F are operators. In this case, since the past motion does not appear, one says 
there is no memory. On the other hand, the autocorrelation of the quantum mechanical 
random force [9, 39] becomes 




(52) 



mx + (x + V'(x) = F(t). 



(53) 




(54) 
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where f2 T = nkT/h. In the limit h — > this becomes the familiar form of classical Brownian 
motion: 



The notion of a Markovian system combines two aspects: a stochastic equation with no 
memory, as in (53), and a delta- function correlation, as in (55). The quantum system is 
never Markovian, in particular the correlation at zero temperature has no delta-function 

The original classical Brownian motion problem corresponds to a free particle (K = 0), 
spectral distribution independent of frequency u> and kT » hu (absence of quantum ef- 
fects), in which case the position auto-correlation function is readily calculated, leading to an 
exact expression for the mean-square displacement, which for large t reduces to the Einstein 
result (4) for the diffusion constant. However, at low temperatures, non-Markovian quan- 
tum effects become important and cannot be neglected. Anomalous diffusion in quantum 
Brownian motion has also attracted much interest as an explanation of various experiments 



Another application of interest relates to Josephson junctions. At first glance one might 
be puzzled as to how the quantum Langevin equation applies in this case. Actually, although 
we have used the language of particle motion in our formulation of this description, it should 
be clear that the description is more general than the language. Thus, the operator x in 
the quantum Langevin equation (6) can be a generalized displacement operator. By this we 
mean an operator x such that a term V(x, t) = —xf(t), with c-number fit), added to the 
microscopic Hamiltonian of the system, results in an added term f(t) on the right-hand side 
of the equation of motion. One can therefore apply this description to an equation which is 
formally similar to the Langevin equation but in which the physical meaning of x is different. 
One must, however, be cautious to check the above generalized displacement property. It 
turns out that, for Josephson junctions, the phase difference of the superconducting wave 
function across the junction is such a generalized displacement coordinate. This enabled us 
to obtain an expression for the power spectrum of the phase fluctuations (generalizing an 
earlier weak coupling limit result of Josephson) as well the mean square deviation of the 
phase and the power spectrum of the voltage fluctuations [41]. 

Turning to the study of small tunnel junctions, here too the quantum Langevin equa- 
tion proved to be an effective tool in determining the mean-square charge fluctuation on 
the junction by treating the charge fluctuation as a generalized coordinate [42]. The quan- 




classical. 



(55) 



[40]. 
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turn Langevin equation has also been used to study quantum transport for a many-body 
system [43], the advantage being that the separation of frictional (dissipative) and random 
(fluctuating) forces gives a natural separation between the conductivity and the noise. 

In general, the separation of fluctuations (noise) and dissipation in the quantum Langevin 
equation (6), together with the fluctuation-dissipation theorems, enables us to systematically 
analyze noise in many different systems. In particular, using (17), we were able to calculate 
the power spectrum of the coordinate fluctuations in a universal model which we presented 
for the detection of noise in gravitational wave detector systems [44] . 

In addition, we found that quantum effects required modification of the famous Onsager 
(classical) regression hypothesis [18], which states that the regression of fluctuations is gov- 
erned by the macroscopic equations describing the approach to equilibrium. In other words, 
the so-called quantum regression theorem is only correct if one makes various approximations 
[45]. In fact, we claim that the correct generalization of the Onsager regression hypothesis 
is the fluctuation-dissipation theorem of Callen and Welton [18] . 

It is also of interest to note that the quantum Langevin equation has been used in a 
very straightforward manner to analyze the so-called Unruh radiation problem. Its use 
in this context was first considered by Sciama and co-workers [46] but they used various 
approximations which left their work open to criticism. However, we carried out an exact 
analysis of an oscillator (the detector) moving under a constant force with respect to zero- 
temperature vacuum and coupled to a one-dimensional scalar field [47]. We showed that, 
contrary to the conclusions reached by Unruh, this system does not radiate despite the fact 
that it thermalizes at the Unruh temperature. 

More recently, the quantum Langevin equation has been used in the general area of meso- 
scopic systems and quantum information. Its use enabled us to incorporate "entanglement 
at all times" , in contrast to the often used but more approximate master equations. In par- 
ticular, it enabled us to analyze exactly the decoherence of a Schrodinger cat superposition 
of Gaussian states, one striking result being the conclusion that decoherence can occur even 
in the absence of dissipation [48] , a result which does not emerge from master equation calcu- 
lations. It should be noted that one finds many different master equations in the literature. 
However, even for what is referred to as the "exact master equation" (which are exact only 
in the sense that they incorporate time-dependent coefficients) [49], we presented what we 
feel is the most transparant exact solution [50, 51]. Our approach started with the Wigner 
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function equivalent equation since the Wigner function provides the same information as the 
corresponding density matrix while making the calculations simpler and more transparent. 
One striking result we obtained is that the exact master equation problem is equivalent to 
the Langevin equation for the initial value problem, which was much easier to solve. We 
found that serious divergences arose at low temperatures [50] and that, even in the high 
temperature regime [52], problems also exist, notably the fact that the density matrix is not 
necessarily positive. In addition, Karrlein and Grabert [49] showed that there is no unique 
master equation, which is connected with the fact that the Onsager regression hypothesis 
fails in the quantum regime. By contrast, our "entanglement at all times" approach, based 
on the use of the quantum Langevin equation, has none of these problems [53]. 

Since entanglement is the essence of various schemes to build small devices for quantum 
information applications (cryptography, quantum computing, etc.), we used the quantum 
Langevin equation to examine how disentanglement occurs in two-body systems. In particu- 
lar, we found that if the temperature is larger than a critical value that disentanglement can 
occur in the absence of dissipation [54] . However, the effect of temperature is very different 
for disentanglement than for decoherence because the effect is constant for all time whereas 
for decoherence the effect increases with time. 

VII. CONCLUSIONS 

We have seen that radiation occurs when a system interacts with a heat bath. It is, in 
essence, a dissipative term (which is not invariant under t — > —t) which occurs in the equation 
of motion of the system. More generally, it occurs in conjunction with fluctuations and the 
relation between these quantities is expressed in the famous quantum fluctuation-dissipation 
theorem. In the absence of a heat bath, we still have quantum fluctuations. However, 
even in the classical case, we have temperatures fluctuations. It is only in the classical 
case where temperature is neglected do we have zero fluctuations so that the dissipative 
(radiation reaction) term only appears. This is mainly the arena for discussing radiation 
effects in electrodynamics. More generally, as we saw above, there are many systems for 
which fluctuations play a key role in addition to radiation reaction. 

We should remark that we have chosen to write the autocorrelation functions in their sym- 
metric form since we found that this choice led to the most elegant presentation. However, 
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clearly other choices are possible, corresponding to the normal and anti-normal ordering of 
the a + (creation) and a (annihilation) operators, but it is clear that in all cases the physical 
implications are the same, a point emphasized by Milonni in his detailed discussion of QED 
effects [55]. 

In our discussion of radiation effects in electrodynamics, we have used as our starting point 
the universally accepted Hamiltonian of quantum electrodynamics. Treating the blackbody 
radiation field as a heat bath, this enabled us to obtain an exact equation of motion of 
an electron in the form of a quantum Langevin equation, in which both dissipation and 
fluctuations appear. Imposing the requirement that the second law of thermodynamics must 
be satisfied leads to an equation of motion that is second order in time with no runaway 
solutions. Thus our equation, in contrast with the well-known Abraham-Lorentz equation, 
is consistent with physical principles. In particular, causality is preserved at all times and 
the optical theorem is satisfied. 
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